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We study the two dimensional Chern insulator and spin Hall insulator on a non-orientable Rie- 
mann surface, the Mobius strip, where the usual bandstructure topological invariant is not defined. 
We show that while the flow pattern of edge currents can detect the twist of the Mobius strip in the 
case of Chern insulator, it can not do so in spin Hall insulator. 

PACS numbers: 



Band insulator with interesting bandstructure topol- 
ogy, the so called "topological insulator", has attracted 
considerable attention lately [1]. Mathematically these 
free-electron insulators are characterized by topological 
invariants in their bandstructure. Examples include the 
Chern number [2, 3] and Z2 index [4-7] for the Chern and 
Z2 insulator, respectively. Physically a hallmark of these 
insulators is their protected boundary states. In two di- 
mensions these are the chiral edge states of the integer 
quantum Hall effect [8, 9], and the helical edge states of 
the spin Hall insulator [10, 11]. In three dimensions the 
boundary states have a massless Dirac fermion dispersion 
relation[12, 13]. 

In this paper we ask a simple question: in two di- 
mensions can we put topological insulators on a non- 
orientable Riemann surface? and if so, what are their 
signatures? The first question was posed to us by Prof. 
X. Sun of the Fudan University. Clearly the usual band- 
structure topological invariants can not be defined in this 
situation; therefore naturally one would examine the edge 
state structure. In the following we study two dimen- 



sional Chern and Z2 insulator on the Mobius strip. Since 
the Mobius strip has only one edge, it is far from clear 
how does the flow pattern of the current look like. 

First we start with the Chern insulator. As pointed 
out by Haldane[14], the necessary ingredient of Chern 
insulator is time-reversal symmetry breaking rather than 
net magnetic flux. Specifically we consider the following 
two-band model on square lattice 



n (k) = d(k) . r. 



(1) 



Here r are Pauli matrices which acts on the (two) orbital 
degrees of freedom, and 



d(k) = (sin/ca^, sin ky^ 1 — coskx — cos ky). 



(2) 



It is straightforward to check that the Chern number (or 
the TKKN index) associated with the valence band of 
this model is 1. In order to implement this model on the 
Mobius strip we first need to Fourier transform the above 
model to real space: 
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In the above (i, j) are the integer coordinates of the sites 
of a square lattice, and ?/; is a two-component fermion 
field associated with the two orbitals in question. Be- 
cause we shall study Eq. (3) on the Mobius strip it is 
essential to define how do the orbitals couple to the local 
orientation. A convenient definition is to let the pseu- 
dospin corresponding to the orbital degrees of freedom 
couple to curvature in the same way the real spin in Dirac 
theory does[15]. This amounts to replacing the r matri- 
ces in Eq. (3) by the following position-dependent Pauli 
matrices, i.e., 



r. 



(4) 



Here /i = 1,2,3 and fiij^^ are unit vectors defining a 
local frame when the relevant surface is embedded in the 
three dimensional Euclidean space (ns is the local surface 
normal). 

Let us warm up by studying the surface of a cylinder. 
We build a coordinate system as follows: 

x{u,v) = cosu^ y{u^v) = sinu^ z{u^v) = v. (5) 

where < u < 27t and — 1 < '^ < 1. Thus we have 
a cylinder of height 2 and radius 1. The local frame is 





FIG. 1: (color on-line) The edge current associated with the 
Chern insulator described in the text. There are 40/12 sites 
in the circumferential and height directions respectively. The 
direction of the current is illustrated by the red arrows. The 
length of arrows is proportional to the magnitude of the cur- 
rent. To make the plot easier to understand we have omitted 
the weak current away from the edge. 



FIG. 2: (color on-line) The edge current associated with the 
Chern insulator on a Mobius strip. There are 40 x 12 sites, 
and all bonds across a line segment are removed. Again, to 
make the plot easier to understand we have omitted the weak 
current away from the edge. 



defined by 



ni{u^v) = dyf{u,v) 
h2{u,v) = dur{u,v) 
hs{u,v) = hi{u,v) X h2{u,v), 



(6) 



where f{u^v) = {cosu^sinu^v). Set u^v to a. set of dis- 
crete values corresponding to the square lattice we sub- 
stitute Eq. (6) into Eq. (4) then Eq. (3). We diagonalize 
the resulting Hamiltonian numerically and compute the 
expectation value of the current operator. 



ji,2 =mjrij.i^2^Pi,j' 



(7) 



In Fig. (1) we plot the result. As expected counter- 
propagating chiral edge currents are found on the two 
opposite edges. 

Now we are ready for the Mobius strip. The coordinate 
system we use is 



c{u,v) = (^ 
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In Fig. (2) we first present the result when all bonds 
across a line segment are removed. We note that a pair 
of co-propagating edge currents are localized on the cut. 
Their presence is due to the orientation flip across the 
cut. As a result the chirality of the Chern insulator re- 
verses across the cut. This is similar to the edge current 
produced at the location of magnetic field reversal in the 




FIG. 3: (color on-line) The edge current associated with the 
Chern insulator on a Mobius strip. There are 40 x 12 sites. 
To make the plot easier to understand we have omitted the 
weak current away from the edge. 



quantum Hall effect. Because the edge currents at the cut 
are co-propagating^ seal the cut has no effect on them. 
The current pattern after the cut is sealed is shown in 
Fig. (3). It is worthy to note that while geometrically 
there is no singularity on a Mobius strip, in defining the 
chirality of the Chern insulator it is necessary to choose a 
cut across which the direction reverses. Thus the current 
pattern of the Chern insulator does detect the twist of 
the Mobius strip. 

Next we study the Z2 insulator on the Mobius strip. 
The momentum space Hamiltonian in M? is given by 

n (k) = di(k)ri (g) /2 + d2{k)T2 ^ h ^ ^3(k)r3 (g) ^3, (9) 

where I2 is the 2x2 identity matrix. The corresponding 



real space version is 






(10) 
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FIG. 4: (color on-line) The edge current associated with the 
Z2 insulator on a Mobius strip. The blue arrows indicate the 
direction of local Sz axis. The red arrows at the top/bottom of 
the blue ones illustrate the current associated with Sz — il, 
respectively. There are 40 x 12 sites, and to make the plot 
easier to understand we have omitted the weak current away 
from the edge. 




FIG. 5: (color on-line) The edge current associated with the 
Z2 insulator on a Mobius strip. Bonds across a line segment 
are removed. The meaning of blue and red arrows are the 
same as in Fig. (4). There are 40 x 12 sites, and to make the 
plot easier to understand we have omitted the weak current 
away from the edge. The inset zooms in at the currents near 
the cut. Note for each spin component there is a pair of 
counter-propagating current. 



In defining Eq. (10) on the Mobius strip we define the lo- 
cal spin Pauli matrices, cr^^j;^, in the same way in Eq. (4). 
We also start by studying the edge current pattern on the 
cylinder. As shown in Fig. (4) there are a pair of time- 
reversal conjugate, counter-propagating edge currents at 
each edge. Here blue arrows indicate the direction of lo- 
cal 5'^ axis. The red arrows at the top/bottom of the 
blue ones illustrate the current associated with Sz = ±1, 
respectively. 

At last we study the Z2 insulator on the Mobius strip. 
Again, we begin by removing the bonds across a line 
segment. The associated current pattern is shown in 
Fig. (5). The meaning of blue and red arrows are the 
same as in Fig. (4). The inset zooms in at the currents 
near the cut. It is important to note that for each spin 
component there is a pair of counter-propagating current. 
As the result, when the cut is sealed they are allowed to 
back scatter against each other and hence gaps out the 
associated edge modes. The result with the cut sealed is 
shown in Fig. (6). As expected the edge currents associ- 
ated with the cut are completely removed. 

In summary we have studied the edge current distribu- 



tion of the Chern and Z2 insulators on the Mobius strip. 
The current pattern of the Chern insulator clearly de- 
tects the twist, while that of the Z2 insulator does not. 
This reveals an interesting interaction between the topol- 
ogy of the electronic structure and that of the substrate. 
The fact that the Z2 topological insulator can be seam- 
lessly put on a Z2 fiber bundle (the Mobius strip is a Z2 
fiber bundle over a circle) is particularly interesting. The 
mathematical meting of this needs to be clarified in the 
future. 
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FIG. 6: (color on-line) The edge current associated with the 
Z2 insulator on a Mobius strip. The meaning of blue and red 
arrows are the same as in Fig. (4). There are 40 x 12 sites, 
and too make the plot easier to understand we have omitted 
the weak current away from the edge. 
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